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ABSTRACT 



We describe the quantum mechanical scattering of slowly moving maximally 
charged black holes. Our technique is to develop a canonical quantization proce- 
dure on the parameter space of possible static classical solutions. With this, we 
compute the capture cross sections for the scattering of two black holes. Finally, 
we discuss how quantization on this parameter space relates to quantization of the 
degrees of freedom of the gravitational field. 



1. Introduction 



Here we study quantum mechanical interactions of charged black holes. In 
classical general relativity there exist exact static solutions for N maximally charged 
black holes; the black holes can be placed anywhere, and will remain at rest. 
This suggests that for slowly moving maximally charged black holes, the spatial 
geometry at any time will be well approximated by the static solution for the 
black hole configuration at that time. The classical solutions were worked out [2] 
following this type of adiabatic approach. Now, the parameter or moduli space of 
possible static solutions is a 3N dimensional manifold, consisting of the positions 
of the N black holes. Therefore, for the slowly moving black holes, a path is traced 
out in moduli space as the 3- geometry evolves to the 4-dimensional spacetime. 

In studying the classical solutions it was found that the motion of the black 
holes was governed by an effective Hamiltonian for N point particles. The approach 
here is to use this Hamiltonian to evolve a Schroedinger wave function, for the case 
of two black holes. The wave function is a function of the positions of the two black 
holes, that is, the configuration space of the wavefunction is the parameter space 
of possible classical static solutions. Hence the degree of freedom being quantized 
is the proper distance between the two black holes. 

One can imagine that if the full theory for quantum gravity were known, one 
could compute the motion of maximally charged black holes and then consider 
the slow motion, low energy limit. Here, we can hope that we are computing 
an approximation to the actual motion. Of course, since we do not know the 
theory of quantum gravity we have no way of knowing how good an approximation 
this is. On the other hand, we do know of some ways this approximation may 
fail. In particular, the classical metrics were assumed to have a particular form, 
corresponding to slow motion. However, in the full quantum theoretic description 
of the system, there will surely be excitation of degrees of freedom which are not 
included in this model, i.e., the configuration space for the true wave function must 
include other degrees of freedom than those in the moduli space considered here. 
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In our approximation we will never see these. This is similar to mini-superspace 
models, and also quantized non-linear sigma models (when viewed as models of 
approximately constrained systems), which suffer the same shortcomingQthere is 
no way excite modes which are not in the approximation. [1] 

Perhaps the most interesting issue is the meaning of the wavefunction for the 
geometry of the spacetime. In the classical case the metric can be found by solving 
constraint equations, with the black hole positions and momenta as sources, after 
the motion of the black holes is known. When the sources are quantized, one can 
no longer speak of a position and velocity of the source, and instead must speak in 
terms of probabilities. The question, then, is what does it mean for the metric field 
equations to have sources which are probabilities? Our preferred interpretation, 
but by no means the only one, is to consider the probabilities as probabilities for 
finding different spacetimes. We discuss these issues in the final section. 

2. The System for Scattering of Maximally Charged Black Holes 

We now review the work of Ferrell and Eardley. [2] This section will conclude 
with the classical Hamiltonian for black holes. In the next section we will discuss 
the quantization procedure, e.g., the Hilbert space of states, and the quantum 
mechanical operator which corresponds to the classical Hamiltonian. 

The sources in our study will be maximally charged black holes. A maximally 
charged black hole is an electrically charged (Reisner-Nordstrom) black hole with 
charge Q = \J~GjM. Such a black hole is maximally charged, since if Q is increased 
but M is held fixed, the event horizon which is associated with the black hole 
disappears and a naked singularity arises. 

The theory for charged matter in curved space is assumed to be the coupled 
Einstein- Maxwell field theory. The action for this theory is [3]: 

S = ^J ~lj F ab F ab + J A a f + J mi dn (2.1) 

In this expression, is the scalar curvature of the spacetime, A a is the electro- 
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magnetic four- potential, j a is the electromagnetic current and F a \, = d a A^ — d\>A a 
is the electromagnetic field strength. The last term is the action for a collection 
of free particles, and doesn't make sense for black holes. In deriving the classical 
effective action, Eardley and Ferrell replace the singularities with charged dust, 
with charge density equal to the mass density. So in performing the manipulations 
the sources are smooth. At the end of the derivation, the limit is taken in which 
the charge density goes to a sum of point singularities. This regularization scheme 
works because the dust is maximally charged, and hence in a static configuration. 
Dust with charge less than the mass would tend to collapse. For this paper we will 
not make the distinction between black holes and the smoothed distribution used 
to represent them. G is Newton's gravitational constant. In this paper we will 
work in units where G = 1 and c = 1. In these units k = ^p/ anc ^ where mp/ anc fc 
is the Planck mass. 

In non-relativistic theories it is clear that there exists a static configuration of 
maximally charged black holes, since the condition for maximal charge is just the 
condition that the static gravitational attraction exactly cancel the electrostatic 
repulsion. Remarkably, in the full Einstein- Maxwell theory for gravitation coupled 
to electricity and magnetism there also can be static configurations of maximally 
charged black holes. 

The metric that describes a system of N static maximally charged black holes 
(first discovered by Majumdar[4] and Papapetrou[5]) is: 

ds 2 = -^- 2 dt 2 + vHijdxidx*, (2.2) 

where V 2 ^ = — 47r£ a m a 5(x — x a ) with the black holes at the points x a and the 
boundary condition ^ = 1 at infinity. Here V 2 is the Laplacian on the x regarded as 
coordinates on flat i? 3 . In this coordinate system the hypersurfaces t = constant 
intersect the singularities at the points x = x a , and the event horizons are the 
same points (i.e., the event horizons are represented as points, but have surface 
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area Anm a .) The Maxwell 4-potential is 

A = -(1 - 1/V)dt. (2.3) 

The space of possible Majumdar-Papapetrou metrics is the 3N dimensional 
configuration space of the positions of the black holes-once the x a are known, \1> 
is known and therefore the metric and scalar potential are also known. If the 
black holes are moving slowly then the solution will trace out a trajectory in 
the space of static solutions. [2,6] These quasi-static solutions exist because in the 
slow motion approximation radiation can be neglected: [3] The static forces cancel 
exactly, and the only relevant forces are the velocity dependent magnetic force 
and the gravitational analog, the gravito-magnetic force. As long as the velocities 
remain small, there will be very little radiation. In this approximation then, the 
fields are fully fixed by the positions and velocities of the matter sources. (A 
similar example of this approximation comes from electrodynamics in flat space. 
If the position of a charge q at ro is changing slowly, with v = dro/dt, then the 
scalar potential is approximately (f>(r,t) ~ q/\r — ro(t)\, and the vector potential 
is approximately A m qv/\r — ro(t)\. The space of field configurations is just the 
configuration space for the particle, since a field configuration is determined by the 
position and velocity of the particle.) 

To describe slowly moving black holes, one looks for a metric and four potential 
of the form 

ds 2 = ~^dt 2 + 2N i dx i dt + ^/ 2 5 ij dx i dx j (2.4) 

A = -(1 - -)dt + Aidx 1 (2.5) 

with ^ as before, but now the x a will be functions of time. Both iV, and A{ 
will depend on the velocities of the black holes. However, since we are assuming 
velocities are small we can truncate the field equations for N{ and A{ to first order 
in velocity. 
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Using the Einstein constraint equations, the fields N{ and Ai can be solved for 
in terms of the source positions and velocities. These expressions are substituted 
into the Einstein-Maxwell action (2.1). The resulting effective Lagrangian depends 
only on the positions and velocities of the sources. The black hole limit is shown to 
be well defined. This is the Lagrangian for the interaction of the black holes. For 
details see [2]. Finally, one finds that the Hamiltonian for two black holes, with 
masses mi and TO2, is 



where P is the momentum of the center of mass, M = mi + nri2, r = xi — x 2 
is the relative coordinate, p is the momentum conjugate to that coordinate, and 
fx = 17111712 /M is the reduced mass of the system. In the first term, a flat metric is 
implied. The metric in the second term is 



where a — 1 — 2fi/M. (Recall that R is the position of the center of mass, and r is 
the relative position, hence these are coordinates on the moduli space of solutions, 
not on spacetime.) 

The evolution of the center of mass coordinate is just free particle motion. The 
evolution for the relative coordinate can have one of two behaviors. If the two black 
holes start at infinite separation, then when they interact they can either scatter 
back out to infinity, or they can evolve toward zero separation, depending on their 
angular momentum. Although in the r — > limit the slow motion approximation 
probably breaks down, it is reasonable to guess that in this case the two black 
holes coalesce into one black hole. [2,3] We will call the two possible classes of orbits 
scattering orbits or coalescence orbits, respectively. This completes the review of 
the classical behavior. 




(2.6) 



9 <ab = l(r)S a b 



(2.7) 



7(0 = 1 + — + — + — 



3M 3M 2 3aM 3 



(2.8) 
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3. Solutions of the Schroedinger equation for Charged Black Holes 



Now we will view the hamiltonian EL<i-body, f° r the system of two maximally 
charged slowly moving black holes, as defining the evolution of a quantum mechan- 
ical system. Hence the system is described by a wave function. As in the classical 
case, the center of mass degrees of freedom X c separate, so that the total wave 
function can be taken to be of the form ^ total = ^( x ) t)exp i(—E c t/h + P c • X c ), 
where E c and P c are the center of mass energy and momentum, and x are the 
relative coordinates. Suppose that at early times a wave packet is given which has 
support at large relative coordinate r, so the two black holes are far apart. As 
time evolves the black holes approach each other, that is the wave packet evolves 
towards smaller r. Part of the packet will be scattered and part will be absorbed; 
classically of course the black holes either scatter or coalsce. Here we will compute 
the absorption coefficient as a function of the angular momentum and energy of the 
wave. The Schroedinger equation for the wave function of the relative coordinates 
is 

"7*(x,f) = - (E c - ^3- - M)* , (3.1) 



where 



and 



h 1 

H = -—g ab V a V b +h\n 
2/i 



/ \ / j 2 2jo2\ , , n 1 3M 3M 2 aM 3 
9ab = l{ r ) [dr +r all ) , and 7(r) = 1 + + H — ^— . 

1Z is the scalar curvature of the three-metric g a b, and a — 1 — 2/i/M. The Hilbert 
space of states can be taken to be square integrable functions. We will see that the 
energy eigenfunctions are square integrable as r — > 0. As r — > oo, the eigenfuntions 
become plane waves, so just as usual, one would need to form wave packets in the 
free particle regime. When checking that H is hermitian, one finds that boundary 
term contributions of probability flux cancel, between large r and the horizon. This 
is similiar to the cancellation for plane waves in a box. 



7 



We take the eigenfunctions to be of the form ^ q i m = ip q i(r)Yi m (Q), so that 
Hty = implies taht the radial wave functions satisfy 



il>"(r) + l—ip'(r) - X V -2fj4Ry(rU= -q 2 l{r)i)(r) (3.2) 

r 27 r l 

where 



2 

9 =^ 

It is useful to change to a new "tortoise" coordinate R, which measures actual 
length along the path. This casts the problem into the form of a standard one- 
dimensional quantum mechanical scattering problem, and the potential is better 
behaved. Let R = J y^dr and let x — r^f 1 l 2 i[). Then the Schroedinger equation 
(3.1) becomes 

x,rr+ (q 2 - V) X = (3.3) 



where 

2^3 h(n,r), r -n >2 )+ ^J- 1 

For ji = we have 



V = ^3 M^r ),r " r 7 ' 2 ) + + 2/i£tt. (3.4) 



3 1 (r/M) 2 /(/ + !) r/M 

2M 2 (l + r/M) 5 M 2 (1 + r/M) 3 l ' J 

Here r is an implicit function of R. The new coordinate R ranges from —00 at 
the horizon to +00; R ~ r + \Mlnj^ for r -> 00 and i? ~ -2 1 /aM(r/M) -1 / 2 as 
r — > 0. For the remainder of the paper we will treat the case £ = 0. We find, valid 
for all /x, 

V, 3 -§ + '^,f or H»M ,3.6) 

^ 24 |! + M±I>, /orfl<< _M. 

Hence the potential falls off rapidly both in the asymptotic free region and as the 
horizon is approached. If we think of the problem as the quantum mechanics of 
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a particle on a curved surface, the three-geometry on which the particle moves 
becomes, as the horizon is approached, flat space minus a three-dimensional wedge 
of solid angle 4n/4 = ir. Hence the eifenfunctions behave like free particles on 
either side of the potential barrier, x — e ±lqR as R — > ±00. We want solutions 
to (3.3) for the eigenf unctions with the boundary conditions that x is purely a 
captured, left-moving wave as R — > —00, so there is no flux out of the horizon. 
Then as R — > 00 the solution will be the sum of an incident wave e~ tqR , normalized 
to unit amplitude, and a scattered wave Se %qR . Using the asymptotic form of the 
potential (3.6), we find 

Xq i~ qR[{-i) l+l hf\qR) + S ql h^\qR)\ , R » M (3.7) 

and 

Xq i ~ C ql qRh?\qR) , i? < -M (3.8) 

where z/ = a/4Z(Z + 1) + 1/4 — 1/2. Hence the fraction of captured flux is \C q i\ 2 
and the scattered flux is \S q i\ 2 = 1 — |C 3 ;| 2 . 

4. General Properties of the Solutions 

First let us consider the motion of wavepackets, so each eigenmode evolves like 

exp(— ihq 2 t/2fi). 

At early tines, let a wave packet start in the asymptotic free region, R — > 00, t — > 
—00. Then by looking at the point of stationary phase, one sees that initially the 
center of the packet moves inward according to R ~ —Vcot, where Vqq = hq/fi 
is the nonrelativistic (relative coordinate) velocity at infinity, associated with the 
momentum fiq. At large times the wave packet has split into two pieces. The center 
of one part continues towards R — > — 00 as R ~ — 2v OQ t; this is the captured flux, 
corresponding to classical coalescence of the black holes. Hence it takes an infinite 
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amount of t-coordinate time for coalescence. There is also a scattered part of the 
packet, propagating as R ~ v^t. Hence the motion of the center of the packet is 
like the classical solution. The motion along the classical path depends only on a 
rescaled time, s = Voot, and so the parameter Voo scales out of the problem. Hence 
the classical solutions depend only on the impact parameter b. However, quantum 
mechanically the quantities of physical interest such as the scattering coefficient, 
will depend on t>oo (or equivalently, q ) as well as the analogue of b. Indeed, from 
the eigenfunction equation (3.3) we see that the eifenfunctions depend on the two 
dimensionless parameters qM and -jfr, where 



A = V^±S (41) 
M qM y ' ' 

or equivalently on I and qM. The definition of jj comes from equating the quantum 
mechanical angular momentum Tiyjl{l + 1) with the classical angular momentum 
fiVoob. Finally, note that the slow motion approximation t>oo <C 1 is equivalent to 

hq < fi (4.2) 

that is, the (quantum mechanical) particle is nonrelativistic. Although we don't 
have sufficient control of the solutions for the eigenfunctions, to actually match the 
two asymptotic regions, we can determine the scaling with qM of the coefficients 
C q i and S q i in the low energy limit. Then in the next section we will use different 
techniques to get approximate forms for C and S in essentially all qM and / 
regimes. Now, suppose qM <C 1 . Then for \R\ < 1/q, the solutions to (3.3) are 
(approximately) independent of q. The solution in the region M < R < 1/q inherits 
an overall scaling dependence on qM from the large R solution (3.7). Hence this 
same scaling with qM is passed onto the solution in the region —1/q < R < —M, 
which can then be matched onto the large negative R region. One finds, for qM <C 1 
and all /, 

S ql ~ H) m + 0(qM) 2l+2 



ql 2^i> + i)r(/ + !) iy ; 



*ll / ( 4 - 3 ) 
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where 7/ is an unknown coefficient which is independent of q. For large / then, the 
capture coefficient is 

M 2 - (^) 6 ' +3 ,1 » 1 ,1M - (4.4) 

and 

|C>| 2 oc (gM) 2 ,/or/ = 0. 

(A similiar analysis in the case of a Klein-Gordon field scattering off a Schwarzchild 
background yields \C q i\ 2 oc (gM) a+2 .) Therefore at low energies the only signifi- 
cant capture occurs for I = waves. 

5. Approximate Methods for finding C q i 

We now turn to the calculation of the capture coefficient. When the energy of 
the particle is well over the potential barrier one can use the Born approximation, 
and the WKB approximation when the particle is well under the barrier, for / ^> 1. 
For I = and low energies, one can use another approximation where the potential 
is replaced by a delta-function. For a given I and qM one has to solve a cubic 
equation to decide if the particle is over or under the barrier. In the case /i = 
this simplifies. If I = 0, then the particle is over the barrier if qM > .23. For I > 1, 
the particle is over the barrier - and hence primarily captured- if qM > | y/l{l + 1), 
which is equivalent to -jfa < 2.5. (Recall that in the classical case, the coalescence 
occurs if the impact parameter is less than 2.5M.) 

a) Low Energies 

In the small qM regime one looks for wave functions of the form \S'\~ 1 l 2 e lS . 
In the first WKB or adiabatic approximation, S = — J R (q 2 — V^^dR', for a 
wave propagating in from R = +00. This neglects terms of order if^H? 2 — V|~ 3 / 2 , 
which are small except near the classical turning points R a and Rb where q 2 = V. 
Further, for the WKB approximation to be valid, the width of the potential Rb — R a 
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must be large compared to the wavelength of the incident particle. For I ^ 



and small incident energies, qM < 2/5y/l(l + 1), then R a ~ -21(1 + 1)/R 2 and 
R b ~ Z(Z + l)/^ 2 , so the WKB approximation is valid if a//(/ + 1) > 2. For / = 
WKB is never valid; one can't fit several wavelengths in under the barrier. Then 
using the standard WKB matching formulae one finds the capture coefficent for 



the case when the particle is well under the barrier, and \Jl(l + 1) 3> 2, 



\C q i\ z = 4( 2 ^ + ^) (5-1) 



where 

R b 

l n e = J dR{V -g 2 )i 

Ra 



In the intervals R a < R < -M and M < R < R b , V goes like l/R 2 , and the 
integration can be done exactly. In the remaining region of integration, q 2 can be 
neglected compared to V, and the value is approximately J^ M dRy/V ~ I (which 
actually turns out to be a fair estimate from doing the integral numerically). Using 
the values for R a and R b given above, and keeping terms of leading order in qM/l, 
we find 

\^ql\ ~ ^ ~ { — ) — <1. ^.^J 

The dependence on gM in (5.2) agrees with the expression which was derived 
from scaling arguments (4.1) (Indeed, one can now approximate the unknown co- 
efficient \ki\ 2 ps l 2 e~ 21 /ir 2 2 G .) The important point is that the capture coefficient 
is very small; performing the integral numerically one finds the capture coefficient 
equals zero to 5 significant figures, in ranges where the WKB approximation is 
reliable. 

For I = and energies below the potential barrier we must use a different 
technique. In this case the wavelength of the incident wave is much greater than 
the width of the potential, as long as qM •< 1, that is, M is small compared to the 
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deBroglie wavelength y / 2fiE/h. (In this regime there is no particular additional 
constraint for consistency with the slow velocity requirement (4.2). It is sufficient, 
that fxM/h = iiM/m^pj . is less than or of order one.) This long wavelength 
wave is not sensitive to the details of the potential and the potential in (3.4) can 
be replaced by a delta-function, V — > LVqS(R), with strength fixed by LVq = 
J dRV(R). 

This problem can be solved, and the capture coefficient is 

ic>i 2 = (i + (^fr 1 « a + (^) 2 ) -1 ' 1 = °> « m « 1 ( 5 - 3 ) 

It is simple to estimate LVq to be a few times 10 _2 /M 2 , and numerical integration 
gives LVq = 2 • 1(T 2 /M 2 . As qM -> the capture coefficient goes like (gM) 2 , 
which agrees with the previous result derived from scaling arguments. 

Classically, if the two black holes approach each other with zero angular mo- 
mentum, they always coalesce. In the quantum mechanical system, for incident 
energies of qM = .1, .01 and .001, the / = mode has a capture coefficient of 
.99, .50 and .001 respectively. The behavior interpolates between almost complete 
capture, as one expects for a particle, to complete scattering to the asymptotic free 
region, as is characteristic of a wave. At very long wavelengths the wave barely 
"sees" the black hole. 

In summary, for incident energies under the barrier we can easily compute 
the captured flux from (5.3), and the higher angular momentum waves are almost 
completely scattered. Only the I = mode has a non-negliable capture coefficient. 
This is not too surprising since an angular momentum barrier is hard to tunnel 
through. ((5.3) can also be used for incident energies over the barrier, .23 < qM < 
1, where it gives a capture coefficient of essentially unity. This agrees with the 
results of the high energy approach discussed next.) 

b) High Energies 

For energies high compared to the height of the potential the simplest approx- 
imation is the Born approximation, which appears to be sufficient here. To this 
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end, we rewrite the eigenfunction equation as an integral equation 



Z(R) = Z°(R) + J dR'G(R, R')V(R')Z(R') : (5.4) 

where £°(R) is any solution to the equation with the potential set to zero, and G(R, 
R') is the Green function for that equation. For the no-flux-out-of-the-horizon 
boundary conditions as before one takes £° = e~ iqR and 

iq(R-R') 

G(R, R') = — — — , forR > R' 



-iq{R~R') 

G(R, R') = — , forR < R'. 

Substituting into (5.4), we have £ as the sum of an ingoing and an outgoing wave, 
and hence can deduce the scattering and capture coefficients, S and C. It is simple 
to check that this approximation does not conserve probability. However, if one 
computes S and C to second order, probability is conserved. Furthermore, only C 
gets a correction at this order. S is already correct to second order, so it suffices 
to compute S to first order and derive \C\ 2 from \C\ 2 = 1 — IS*] 2 . 

The range of validity of the Born approximation is 



qM > .23, if 1 = 



qM > (2/5)y/l(l + l) if I > 0. 

In these ranges one has 

|^| 2 = ^2l / dRV(R)e-** R \ 2 (5.5) 

The integral is approximately zero where the Born approximation is valid: the 
width of the potential at half-maximum is about 3M, which means that for qM ^> 
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1, the exponential in the integrand is rapidly oscillating and different contributions 
cancel. We computed the integral numerically and indeed one finds that, for qM > 
10, once the particle is over the barrier, the scattering coefficient is zero to five 
significant figures. 

c) Summary of Scattering Behavior 

For qM > 10 we see particle like behavior. Recall that the condition for the 
wave to be over the barrier, and hence captured, is 



± = vwm < 2 . 5> 

M qM ~ ' 

which is the same as the classical condition for coalescence. For high enough ener- 
gies the discreteness of 1 is no longer important; a sketch of the capture coefficient 
vs. 1 looks very much like the plot of capture vs. a continuous b/M. Thus, for 
high energies, \C\ 2 ~ 1 for b/M < 2.5. For b/M > 3, where we know the WKB 
approximation is valid, \C\ 2 ~ 0. Indeed, \C\ 2 is probably close to zero at an 
impact parameter closer to 2.5 than 3. 

At intermediate energies, .23 < qM < 10, q 2 is still above the potential barrier 
for some /—values. As /, or b, is increased, \C\ 2 decreases from (nearly) unity to 
(nearly) zero. The transition from the Born regime, where \C\ 2 « 1, to the WKB, 
where \C\ 2 ~ 0, is not as sharp as in the high energy case. For qM = .5 this 
transition occurs over a change in b of about 14M, compared to M/2 above. For 
qM < .23, we see transition to wave-like behavior, as previously discussed. The 
wave is almost completely scattered, the captured flux going to zero like (qM) 6l+3 
for large 1 and as (qM) 2 for I = 0. 

Apriori, one might expect more back scattering in the over-the-barrier case 
and more capture in the under-the-barrier case than we have found; except for the 
lowest angular momentum modes, the transition from a capture coefficient of 1 to 
is quite abrupt in b/M. This is due to the "featurelessness" of the black hole 
potential. There is only one length scale, M, which determines both the height 



15 



and the width of the potential. Typically in scattering problems there are two 
independent parameters to vary. 

One can also think of the over-the-barrier condition as fixing q and letting M 
increase, since the analysis is valid in a strong gravity regime. Then, as one expects 
for a black hole, once particle flux is in, it never gets out. Note that because of 
the slow motion approximation (4.2), in the regime qM > 1 one can only consider 
black holes with masses such that 

— + —2 >!» ( 5 - 6 ) 

planck 

in the case when the particle is over the barrier. 



6. Discussion 

At this point the most interesting question is, what has been quantized? When 
solving the classical problem one first solves for the trajectories of the black holes 
which immediately gives the diagonal metric components and the time component 
of the Maxwell potential. Then one finds (at least in principle) the other field 
components by solving given constraint equations, in which the black hole positions 
and velocities appear as source terms. [2] The complete solution can be found in 
this sequence of steps because of the slow motion approximation. 

In the present calculation we have quantized the distance between the two 
black holes. Instead of having a solution r(t) which is this distance as a function 
of time, one now has an amplitude ty(r,t). 

This is fundamentally different from calculations in which "gravitons" are 
quantized as linear perturbations off some fixed classical background (typically 
Minkowski or Robertson Walker). In the charged black hole problem, one of the 
degrees of freedom of the full metric has been quantized, not just a fluctuation. 
To see this, just write the spacetime metric (2.4)in coordinates such that one of 
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the black holes is at the origin and the other is at comoving coordinate z, so that 
one of the metric components is r(t) = a(t)z, the distance between the two black 
holes. 

On the other hand, we have not allowed for quantum fluctuations in " directions 
in the space of metrics" other than this single function. In a fully quantum me- 
chanical theory one must allow all the field components to vary from the classical 
values. The situation is similar to the quantization of a nonlinear sigma model. 
Suppose we have a classical theory in D dimensions (or with D fields) with a strong 
potential that approximately constrains the dynamics to a D - 1 dimensional sub- 
manifold with metric g^. This leads us to model the system by the sigma model 
L = g ab V 'a0V&0, where is constrained to take values on the D - 1 dimensional 
submanifold. Quantizing the sigma model also misses out on quantum fluctua- 
tions in directions off the submanifold. One expects that quantizing the reduced 
system is a good approximation to the full system, if the constraining potential is 
sufficiently strong. In the black hole problem this corresponds to the slow motion 
approximation remaining good. 

An alternative way to say this, is that the approximation made here is the 
same as in "mini superspace" models. In these quantum cosmology models, the 
wave function is taken to depend on only one metric component, namely the scale 
factor. 

We also note that in deriving our black hole Hamiltonian, H e ff, the classical 
equations of motion have been inserted into the action. Now, it was checked 
by Eardly that to derive the classical effective Hamiltonian, one could either work 
exclusively with the equations of motion, or with the action and equations of motion 
for some of the fields, as described earlier. However, since all field configurations 
contribute to a functional integral, one expects that use of the classical equations 
for Ni and A{ also leads to differences from a full quantum theory. Again, this is 
similar to the sigma model approximation discussed above. 

Next, what about the other fields in the problem? We suggest that the prescrip- 
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tion for recovering the metric and electro-magnetic fields from the wave function for 
the black hole coordinates is similar to the prescription for recovering the states 
of Schroedinger's cat in that famous demonstration of quantum mechanics. [11] 
That is, observed macroscopic states occur with probabilities predicted by quan- 
tum mechanics, but are not superpositions of multiple quantum mechanical states. 
Suppose for particular initial conditions of the wave function there is a scattering 
coefficient, p. Then at late times, according to this prescription, one observes fields 
due to two widely separated black holes with masses mi and mi with probability 
p, and the set of fields corresponding to one large black hole with mass M, with 
probability 1 - p. That is, one predicts that certain field configurations occur with 
certain probabilities. 

An alternative prescription would be to couple the classical fields to the ex- 
pectation value of the charge operator. Then, e.g., x would be the sum of two 
Coulombic pieces, with widely separated poles, one with strength m\ + 777.2(1 — p) 
and the other with strength m,2P- That is, the resulting fields are prescribed de- 
terministically, and correspond to sources where "part of the black holes coalesced 
and part scattered", rather than different classical configurations occurring with 
different probabilities. 

Clearly in the gravitational problem one can only argue by analogy and accord- 
ing to what we observe at accessible energies, to choose the theory for specifying 
the classical fields. At this point in our understanding of quantum gravity, we leave 
the choice to the sense and sensibility of the reader. 

Acknowledgements: We would like to thank Doug Eardley and Gary Horowitz for 
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